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Abstract 

We implement the mechanism of non-thermal leptogenesis in the framework of an 
inflationary model based on a supersymmetric (SUSY) Pati-Salam Grand Unified Theory (GUT). 
In particular, we show that inflation is driven by a quartic potential associated with the Higgs 
fields involved in the spontaneous GUT symmetry breaking, in the presence of a non-minimal 
coupling of the inflaton field to gravity. The inflationary model relies on renormalizable 
superpotential terms and does not lead to overproduction of magnetic monopoles. It is 
largely independent of one-loop radiative corrections, and it can be consistent with current 
observational data on the inflationary observables, with the GUT symmetry breaking scale 
assuming its SUSY value. Non-thermal leptogenesis is realized by the out-of-equilibrium 
decay of the two lightest right-handed (RH) neutrinos, which are produced by the inflaton 
decay. Confronting our scenario with the current observational data on light neutrinos, the 
GUT prediction for the heaviest Dirac neutrino mass, the baryon asymmetry of the universe 
and the gravitino limit on the reheating temperature, we constrain the masses of the RH 
neutrinos in the range (10 10 — 10 15 ) GeV and the Dirac neutrino masses of the two first 
generations to values between 0.1 and 20 GeV. 
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1 . Introduction 

One of the most promising and well-motivated mechanisms for the generation of the Baryon 
Asymmetry of the Universe (BAU) is via an initial generation of a lepton asymmetry, which can 
be subsequently converted to BAU through sphaleron effects - see e.g. Ref. (HE). Non-Thermal 
Leptogenesis (nTL) J3J [4j is a variant of this proposal, in which the necessitated departure 
from equilibrium is achieved by construction. Namely, the right-handed (RH) neutrinos, v\, 
whose decay produces the lepton asymmetry are out-of-equilibrium at the onset, since their 
masses are larger than the reheating temperature. Such a set-up can be achieved by the direct 
production of v c t through the inflaton decay, which can also take place out-of-equilibrium. 
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Therefore, such a leptogenesis paradigm largely depends on the inflationary stage, which it 
follows. 

In a recent paper (5]l - for similar attempts, see Ref. (6-8 1 -, we investigate an inflationary 
model where a Standard Model (SM) singlet component of the Higgs fields involved in the 
spontaneous breaking of a supersymmetric (SUSY) Pati-Salam (PS) Grand Unified Theory (GUT) 
can produce inflation of chaotic-type, named non-minimal Higgs Inflation (nMHI), since there 
is a relatively strong non-minimal coupling of the inflaton field to gravity i9Hl2t. This GUT 
provides a natural framework to implement our leptogenesis scenario, since the presence of 
the S(I(2) R gauge symmetry predicts the existence of three 1/ C . In its simplest realization this 
GUT leads to third family Yukawa unification (YU), and does not suffer from the doublet-triplet 
splitting problem since both Higgs doublets are contained in a bidoublet other than the GUT 
scale Higgs fields. Although this GUT is not completely unified - as, e.g., a GUT based on 
the SO(10) gauge symmetry group - it emerges in standard weakly coupled heterotic string 
models (13 1 and in recent D-brane constructions flit . 

The inflationary model relies on renormalizable superpotential terms and does not lead to 
overproduction of magnetic monopoles. It is largely independent of the one-loop radiative 
corrections (15), and it can become consistent with the fitting 1 16 1 of the seven-year data of 
the Wilkinson Microwave Anisotropy Probe Satellite (WMAP7) combined with the baryon-acoustic 
oscillation (BAO) and the measurement of the Hubble constant (Hq). At the same time the GUT 
symmetry breaking scale attains its SUSY value and the \i problem of the Minimal SUSY SM 
(MSSM) is resolved via a Peccei-Quinn (PQ) symmetry, solving also the strong CP problem. 
Inflation can be followed by non-thermal leptogenesis, compatible with the gravitino (G) limit 
tl"7HT9t on the reheating temperature, leading to efficient baryogenesis. In Ref. ((5) we connect 
non-thermal leptogenesis with neutrino data, implementing a two-generation formulation 
of the see-saw i20H22l mechanism and imposing extra restrictions from the data on the 
light neutrino masses and the GUT symmetry on the heaviest Dirac neutrino mass. There 
we assume that the mixing angle between the first and third generation, 13 , vanishes. 
However, the most updated (23 : 24) analyses of the low energy neutrino data suggest that 
non-zero values for 13 are now preferred, while the zero value can be excluded at 8 standard 
deviations. Therefore, a revision of our results, presented in Ref. (5J, is worth pursuing. 

The three-generation implementation of the see-saw mechanism is here adopted, following a 
bottom-up approach, along the lines of Ref. (25-28). In particular, we use as input parameters 
the low energy neutrino observables considering several schemes of neutrino masses. Using 
also the third generation Dirac neutrino mass predicted by the PS GUT, assuming a mild 
hierarchy for the two residual generations and imposing the restriction from BAU, we 
constrain the masses of vf's and the residual neutrino Dirac mass spectrum. Renormalization 
group effects (28: 29 ) are also incorporated in our analysis. 

We present the basic ingredients of our model in Sec. [2] In Sec. [3] we describe the inflationary 
potential and derive the inflationary observables. In Sec. [4] we outline the mechanism of 
non-thermal leptogenesis, while in Sec. [5] we exhibit the relevant imposed constraints and 
restrict the parameters of our model. Our conclusions are summarized in Sec. [6] Throughout 
the text, we use natural units for Planck's and Boltzmann's constants and the speed of light 
(ft = c = fcg = 1); the subscript of type ,x denotes derivation with respect to (w.r.t) the field x 
(e.g., iXX = d 2 /dx 2 ); charge conjugation is denoted by a star and log [In] stands for logarithm 
with basis 10 [e] . 
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2. The Pati-Salam SUSY GUT Model 

In this section, we present the particle content (Sec. |2.1| , the structure of the superpotential 
and the Kahler potential(Sec. |2.2| and describe the SUSY limit (Sec. |2.3| of our model. 



2.1 Particle Content 

We focus on a SUSY PS GUT model described in detail in Ref. (5||30). The representations and 
the transformation properties of the various superfields under Gpg = SlZ(4)c X SIT(2)l X 
S(J(2)r, their decomposition under G$m = Slf(3)c X S(J(2)l x lf(l)y, as well as their extra 
global charges are presented in Table [T] 

The z'th generation (i = 1,2,3) left-handed (LH) quark and lepton superfields, u ia , d la (a = 
1,2,3 is a color index), e; and v, are accommodated in a superfield F,-. The LH antiquark 
and antilepton superfields uf , d c i& , and vf are arranged in another superfield Ff. The 
gauge symmetry Gps can be spontaneously broken down to G$m through v.e.vs which the 
superfields H c and H c acquire in the directions v c H and v c H . The model also contains a gauge 
singlet S, which triggers the breaking of Gpg, as well as an SU(4)c 6-plet G, which splits 
into g c a and g c a under G$m and gives fL3t superheavy masses to d c Ha and d c Ha . In the simplest 
realization of this model 1 13 : 301, the electroweak doublets H u and H rf , which couple to the up 
and down quarks respectively, are exclusively contained in the bidoublet superfield H. 

In addition to Gpg, the model possesses two global (1(1) symmetries, namely a PQ and an R 
symmetry, as well as a discrete Z™ p symmetry ('matter parity') under which F, F c change sign. 
The last symmetry forbids undesirable mixings of F and H and/ or F c and H c and ensures the 
stability of the lightest SUSY particle (LSP). The imposed !i(l) R symmetry, [i(l) R , guarantees 
the linearity of the superpotential w.r.t the singlet S. Finally the U(l) PQ symmetry (J(I)pq, 
assists us to generate the ^-term of the MSSM. The PQ breaking occurs at an intermediate scale 
through the v.e.vs of P, P, and the pterin is generated via a non-renormalizable coupling of 
P and H. Following Ref. l30t , we introduce into the scheme quartic (non-renormalizable) 
superpotential couplings of H c to Ff, which generate intermediate-scale masses for the 
and, thus, masses for the light neutrinos, v,-, via the seesaw mechanism j20H22t . Moreover, 
these couplings allow for the decay of the inflaton into vf, leading to a reheating temperature 
consistent with the G constraint with more or less natural values of the parameters. As shown 
finally in Ref. {30}, the proton turns out to be practically stable in this model. 



2.2 Superpotential and Kahler Potential 

The superpotential W of our model splits into three parts: 

W = W m ssm + W PQ + W HPS , (1) 
which are analyzed in the following. 



Wmssm is the part of W which contains the usual terms - except for the ]i term - of the 
MSSM, supplemented by Yukawa interactions among the left-handed leptons and vf: 



Wmssm = y^HFj = 
y,j (H/eLtf - HjeLtf + H d J eQ la d c j& - H u T eQ ia u c ja ) , with e 



1 
-1 



(2) 



4 



Leptogenesis and Neutrino Masses in an Inflationary SUSY Pati-Salam Model 



Super- 


Represe- 


Trasfor- 


Decompo- 




Global 




fields 


ntations 


MATIONS 


sitions 


Charges 




under G PS 


UNDER G PS 


UNDER G SM 


R 


PQ 2 




Matter Superfields 




(4, 2,1) 


Fi Ul U J 


Q,a(3,2,l/6) 


1 


-1 










Li (1,2, -1/2) 








n 


(4, 1,2) 


U C U R F i 


4,(3, 1,-2/3) 


1 













^(5,1,1/3) 














if (1,1,0) 






















Higgs Superfields 


H c 


(4, 1,2) 




u c Ha (3, 1,-2/3) 








+ 








d e Ha (3, 1,1/3) 














v^j (1,1,0) 














eff (1,1/1) 










(4, 1,2) 


RCU R U C 


«Ha (3/ 1/2/3) 








+ 








^(3,1,-1/3) 














^(1,1,0) 














^(1,1, -1) 








S 


(1, 1,1) 


S 


S(l, 1,0) 


2 





+ 


G 


(6,1,1) 


u c guI 


£(3,1,-1/3) 


2 





+ 








£(3,1,1/3) 








H 


(1, 2,2) 


U L HUl 


H« (1,2, 1/2) 





1 


+ 








H d (l,2,-l/2) 








P 


(1, 1,1) 


P 


P (1,1,0) 


1 


-1 


+ 


P 


(1, 1,1) 


P 


P(1,1,0) 





1 


+ 



Table 1 . The representations, the transformations under Gpg, the decompositions under GgM as well as 
the extra global charges of the superfields of our model. Here U c £ SU(4) C/ U L £ SU(2) L/ (I R £ SU(2) R 
and T, t and * stand for the transpose, the hermitian conjugate and the complex conjugate of a matrix 
respectively. 

Here Q, a = ^w, a d ia j and L, = e ( j are the i-th generation S(I(2) L doublet LH 
quark and lepton superfields respectively. Summation over repeated color and generation 
indices is assumed. Obviously the model predicts YU at Mqut since the fermion masses 
per family originate from a unique term of the PS GUT. It is shown (31 32 1 that exact third 
family YU combined with non-universalities in the gaugino sector and/ or the scalar sector 
can become consistent with a number of phenomenological and cosmological low-energy 
requirements. On the other hand, it is expected on generic grounds that the predictions 
of this simple model for the fermion masses of the two lighter generations are not valid. 
Usually this difficulty can be avoided by introducing d33t an abelian symmetry which 
establishes a hierarchy between the flavor dependent couplings. Alternatively the present 
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model can be augmented 1 34 1 with other Higgs fields so that H u and H d are not exclusively 
contained in fi, but receive subdominant contributions from other representations too. 
As a consequence, a moderate violation of exact YU can be achieved, allowing for an 
acceptable low-energy phenomenology, even with universal boundary conditions for the 
soft SUSY breaking terms. However, we prefer here to work with the simplest version 
of the PS model, using the prediction of the third family YU in order to determine the 



corresponding Dirac neutrino mass - see Sec. 5.1 



WpQ, is the part of W which is relevant for the spontaneous breaking of IT(1)pq and the 
generation of the term of the MSSM. It is given by 

W PQ = A PQ ^ - A„-^-Tr (HeH^e) , (3) 



M s ' 2M S 

where Mg ~ 5 ■ 10 17 GeV is the String scale. The scalar potential, which is generated by the 
first term in the RHS of Eq. j3), after gravity-mediated SUSY breaking, is studied in Ref . ( 30 
[35). For a suitable choice of parameters, the minimum lies at | (P) \ = \ (P) | ~ ^/m 3 / 2 Ms. 
Hence, the PQ symmetry breaking scale is of order ^/m 3 / 2 Ms ~ (10 10 — 10 11 ) GeV. The 
^-term of the MSSM is generated from the second term of the RHS of Eq. J5| as follows: 



A ''2M^ Tr ( H£HTf ) =l iH d J £H u => F ~A,,^-, (4) 



which is of the right magnitude if A fl ~ (0.001 — 0.01). Let us note that V P q has an 
additional local minimum at P = P = 0, which is separated from the global PQ minimum 
by a sizable potential barrier, thus preventing transitions from the trivial to the PQ vacuum. 
Since this situation persists at all cosmic temperatures after reheating, we are obliged to 
assume that, after the termination of nMHI, the system emerges with the appropriate 
combination of initial conditions so that it is led 1 36 1 in the PQ vacuum. 

W HP g, is the part of W which is relevant for nMHI, the spontaneous breaking of Gpg and 
the generation of intermediate Majorana [superheavy] masses for v\ [d c H and d c H ]. It takes 
the form 

W hps = AS [H C H C - Mp s J + \ H H cJ GeH c + \ n H c GeH cJ + k iv c M ^ > (5) 



where Mps is a superheavy mass scale related to Mqut _ see Sec. 3.2 - and G is the dual 



tensor of G. The parameters A and Mpg can be made positive by field redefinitions. 

According to the general recipe fTTl[T2t . the implementation of nMHI within SUGRA requires 
the adoption of a Kahler potential, K, of the following type 

n 2, (\ HCf H c H C H C+ Tr (G + G) |S| 2 , |S| 4 k H ,- CTrC . A 

K = -3m P ln 1 - — — - i-L +k s ^\ + -*L {H H L + h.c.) , 

V 3mp 3mp 6mp 3mp 3mp 2mp J 

(6) 

where nip = 2.44 ■ 10 18 GeV is the reduced Planck scale and the complex scalar components 
of the superfields H c , H c , G and S are denoted by the same symbol. The coefficients kg and k H 
are taken real. From Eq. |6j we can infer that we adopt the standard quadratic non-minimal 
coupling for Higgs-inflaton, which respects the gauge and global symmetries of the model. 
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This non-minimal coupling of the Higgs fields to gravity is transparent in the Jordan frame. 
We also added the fifth term in the RHS of Eq. |6) in order to cure the tachyonic mass problem 
encountered in similar models iTMl2t - see Sec. 3.1 In terms of the components of the various 
fields, K in Eq. |6| reads 



K ■■ 



-3mp In 



3m 2 



+ S 3m| + 2m 2 



v H v H + e H e H + u H a H 



H d c H +Kc) 



(7a) 



with 

f = v H,Vwe c H ,? H ,uj 1 ,a c H ,d c H ,d c H ,g c ,g c and S 
and summation over the repeated Greek indices is implied. 



(7b) 



2.3 The SUSY Limit 

In the limit where m P tends to infinity we can obtain the SUSY limit of the SUGRA potential. 
Assuming that the SM non-singlet components vanish, the F-term potential in this limit, V P , 
turns out to be 

V ¥ = A 2 \v c H v c H - Mp s | 2 + A 2 |S| 2 (Jt/^l 2 + \v c H \ 2 ^ , (8a) 
while the D-term potential is 

Vd=%(\v c h\ 2 -\* c h\ 2 ) 2 ■ (8b) 

Restricting ourselves to the D-flat direction \v c H \ = \v c H \, we find from Vp that the SUSY 
vacuum lies at 

<S)~0 and \{v c H )\ = \{v c H )\ =M PS . (9) 

Therefore, Whps leads to spontaneous breaking of Gps- As we shall see in Sec. [3} the same 
superpotential, W HP g, gives rise to a stage of nMHI . Indeed, along the D-flat direction 
\ V H = \ v h\ anc ^ S = 0, VgusY tends to a quartic potential, which can be employed in 

conjunction with K in Eq. |6f for the realization of nMHI along the lines of Ref. fl2t . 

It should be mentioned that soft SUSY breaking and instanton effects explicitly break U(1)r x 
LJ(1)pq to Z2 x Zg. The latter symmetry is spontaneously broken by (P) and (P) . This would 
lead to a domain wall problem if the PQ transition took place after nMHI. However, as we 
already mentioned above, (J(I)pq is assumed already broken before or during nMHI. The 
final unbroken symmetry of the model is G^m x Z™ p . 



3. The Inflationary Scenario 

Next we outline the salient features of our inflationary scenario (Sec. \3.1\ and calculate a 
number of observable quantities in Sec.|3.2| 



3.1 Structure of the Inflationary Potential 

At tree-level the Einstein Frame (EF) SUGRA potential, Vm, is given by fiTt 

V m = - 3^1) + \i £ D„D fl , (10a) 
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where g is the unified gauge coupling constant and the summation is applied over the 21 
generators T„ of the PS gauge group - see Ref. 1 5 ). Also, we have 

Kofi = Vff' KhK ^ = 4' F « = W HPV« + V w hps/»4 ^d D a = <p K (T„)« (10b) 
The (p a 's are given in Eq. |7b| . If we parameterize the SM singlet components of H c and H c by 

v c H = he' e cos 6 V I Vl and v c H = he iB sin V / \fl, (11) 
we can easily deduce that a D-flat direction occurs at 

6 = 9 = 0, 6 V = tt/4 and e c H = <? H = u c H = u c H = d c H = d c H = g c = f = 0. (12) 



Along this direction, the D-terms in Eq. 1 10a < - and, also, Vq in Eq. |8bJ - vanish, and so Vhi 
takes the form 



with 



/ = 1 + c n x\, m PS 



Mps 
nip 



\ 2 (x 2 h -4m FS ) 2 
16/ 2 

h 1 k H 

■ x >' = and c k = --z + -T- 

fflp o 4 



(13) 



(14) 



From Eq. 1 13 , we can verify that for c-ji ^> 1 and m P g <C 1, Vhi takes a form suitable 
for the realization of nMHI, since it develops a plateau - see also Sec. 3.2 The (almost) 
constant potential energy density V H io and the corresponding Hubble parameter H H i (along 
the trajectory in Eq. (T2)) are given by 



A 2 /x 4 A 2 m p 

Vhtd — 



V, 



1/2 



16/ 2 



andHH I = ^^^ 

16c^ V3m P 4V3. 



(15) 



We next proceed to check the stability of the trajectory in Eq. (T2) w.r.t the fluctuations of the 
various fields. To this end, we expand them in real and imaginary parts as follows 



xi + tx 2 v x\ + ix 2 



where X = e c H , u c H , d c H , g c and x = e, u, d, g . 



(16) 



Notice that the field S can be rotated to the real axis via a suitable R transformation. Along the 
trajectory in Eq. (T2| we find 



k *b = dia e 



M K 1 1 

3+6-3 times 



with Mj; = 



K K 
K K 



, k = 3^x1 and k = f + K. (17) 



To canonically normalize the fields v c H and v c H , we first diagonalize the matrix Mj(. This can 
be achieved via a similarity transformation involving an orthogonal matrix Uk as follows: 



U K M K ul = diag (/,/) , where / = / + bc\x\ and U K 



V2 



1 1 
-11 



(18) 
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Utilizing U^, the kinetic terms of the various fileds can be brought into the following form 

= ^2 ( h2 + \ hW +) + Yf {¥ l + + 2j* K *« = W + (19) 

where 0± = (9 ± 9) / sjl, 

Xa — Xi,X2,Xi,X2,S and ip a = 9 + ,9-,9 v ,Xa an d the dot denotes 
derivation w.r.t the cosmic time, t. In the last line, we introduce the EF canonicaiiy normalized 
fields, h and tp, which can be obtained as follows - cf. Ref. dSl fTTl [1211371 : 

dh f y 2 s/lf v7 V 4 > V? 

Taking into account the ap prox imate expressions for h, } and the slow-roll parameters e, rj, 
which are displayed in Sec. |3.2| we can verify that, during a stage of slow-roll inflation, 9 + ~ 
Jh9+ 1 y/2 since Jh ~ V6m F , 0_ ~ h9_ / y/Zf and 9 V ~ h9 v / y/f since h/y/f ~m P / y/c^. On 
the other hand, we can show that Xa — Xa/ \ff > since the quantity //2/ 3 ' 2 ^ a/ involved in 
relating Xa to Xa, turns out to be negligibly small compared with Xa- Indeed, the Xa's acquire 
effective masses mo S> Hm - see below - and therefore enter a phase of oscillations about 
Xa = with decreasing amplitude. Neglecting the oscillatory part of the relevant solutions, 
we find 

X ~ XaoVfe~™ /3 ^d ~ -2x^]Hnm^ n \ (21) 

where £„o represents the initial amplitude of the oscillations, rj Xli = m~ / 3Hhi and we assume 

Xa(t = 0) = 0. Taking into account the approximate expressions for h and the slow-roll 
parameter e in Sec. |3.2| we find 



- f/2f /2 x. = [c n eHyml) % R « (22) 
Having defined the canonicaiiy normalized scalar fields, we can proceed in investiga ting the 



stability of the inflationary trajectory of Eq. 12 <. To this end, we expand Vhi m Eq. JlOal to 



quadratic order in the fluctuations around the direction of Eq. (12) , as described in detail in 
Ref. (5). In Table[2]we list the eigenvalues of the mass-squared matrices 



with ip x = 9+,8-,9 v ,x\,X2,x\,X2 and S (23) 

Eq.p) 



involved in the expansion of Vhi- We arrange our findings into three groups: the SM singlet 
sector, S — v c H — v c H , the sector with the u c H , u c H and the e c H , e c H fields which are related with 
the broken generators of Gpg and the sector with the d c H ,d c H and the g c ,g c fields. Upon 
diagonalization of the relevant matrices we obtain the following mass eigenstates: 

x\± = —= [x\ ± Xi) and £2^ = — p fe T £2) with z = u, e, d and g. (24) 
v2 V2 

As we observe from the relevant eigenvalues, no instability - as the one found in Ref. (37 1 - 
arises in the spectrum. In particular, it is evident that k$ > 1 assists us to achieve m~ > 
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Table 2. The scalar mass spectrum of our model along the inflationary trajectory of Eq. \12f . To avoid 
very lengthy formulas we neglect terms proportional to nip S and we assume Ah — f° r tne derivation 
of the masses of the scalars in the superfields d c H and d c H . 

- in accordance with the results of Ref. (12 1. Moreover, the D-term contributions to ml and 

ml_ - proportional to the gauge coupling constant g ~ 0.7 - ensure the positivity of these 

masses squared. Finally the masses that the scalars d\ 2 acquire from the second and third 
term of the RHS of Eq. |5| lead to the positivity of ml for A H of order unity. We have also 

numerically verified that the masses of the various scalars remain greater than the Hubble 
parameter during the last 50 — 60 e-foldings of nMHI, and so any inflationary perturbations 
of the fields other than the inflaton are safely eliminated. 

The 8 Goldstone bosons, associated with the modes x\ + and xi- with x = u a and e, are not 
exactly massless since VfD ^ - contrary to the situation of Ref. 1 30 1 where the direction with 
non vanishing {v c H ) minimizes the potential. These masses turn out to be m X Q = Ampx h /2f. 
On the contrary, the angular parametrization in Eq. (TT) assists us to isolate the massless mode 
9- , in agreement with the analysis of Ref. 1 1 1 1 . Employing the well-known Coleman- Weinberg 
formula <15> , we can compute the one-loop radiative corrections to the potential in our model. 
However, these have no significant effect on the inflationary dynamics and predictions, since 
the slope of the inflationary path is generated at the classical level - see the expressions for e 
and ff below. 

3.2 The Inflationary Observables 



Based on the potential of Eq. 1 13 i and keeping in mind that the EF canonically inflaton h is 
related to h via Eq. (20), we can proceed to the analysis of nMHI in the EF, employing the 
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standard slow-roll approximation. Namely, a stage of slow-roll nMHI is determined by the 
condition - see e.g. Ref. (38: 39): 



max{e(7i), \tf(h)\} < 1, 



where 



and 



m l ( V mji \ _ m p ( v Hi,h \ ^ ^ imp 



Vm J 2 / 2 1 V m I 3c\\fi 



(25a) 



~_ 2 V m,hh _ m F ( v ni,hh V mih J ih \ 4:f mp 



are the slow-roll parameters and fo — f ((h) — 2Mpg) = 1 + 4cnmp 5 - see Sec. 
employ Eq. Jl5j and the following approximate relations: 



4.1 



Here we 



/ ~ v#£, V mJl ~ ^%hml and V mJlh ~ - 12 ^/ /np. (26) 

The numerical computation reveals that nMHI terminates due to the violation of the e criterion 
at a value of h equal to hf, which is calculated to be 

e(h f ) =l=>hf= (4/3) 1/4 m P vToT^. (27) 



The number of e-foldings, N*, that the scale fc* = 0.002/Mpc suffers during nMHI can be 
calculated through the relation: 

a. = 4 /Nsjk * (28) 



where is the value of h [h] when fc* crosses the inflationary horizon. Given that h f <§C h*, 

we can write h* as a function of as follows 

3c K hi - h\ 
4/ m p 



N* ~ -t^^t- 1 => ^* = 2m P JN,f /3c n - (29) 



The power spectrum of the curvature perturbations generated by h at the pivot scale fc* is 
estimated as follows 

1 y m(i) 3/2 ^ Afc2 ^ AN, . (30) 



2V3mnp |t? m ^(/j*)l 16V2nf m F 12y/2nc n 



Since the scalars listed in Tableware massive enough during nMHI, A7J can be identified with 
its central observational value - see Sec. [s]- with almost constant N*. The resulting relation 
reveals that A is to be proportional to c^. Indeed we find 



A ~ 8.4 ■ KT 4 ttck/N* => c n ~ 20925A for N* ~ 55. 



(31) 
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The (scalar) spectral index n s , its running a s , and the scalar-to-tensor ratio r can be estimated 
through the relations: 

n s = 1 - 6e* + 2tf* ~ 1 - 2/N*, (32a) 



and 



«s = 3 (4^ - (n s - 1) Z J - 2^ ~ -2£* ~ -2/N^ (32b) 

r = 16e* ~ 12/N^, (32c) 
where £ = m p^Hi,ft^ HI ^/^Hi = m F \/2erfjj/ } + 2tfe. The variables with subscript * are 
evaluated at h = h* and Eqs. |25a) and \25b\ have been employed. 

4. Non-Thermal Leptogenesis 

In this section, we specify how the SUSY inflationary scenario makes a transition to the 
radi ation dominated era (Sec. |4.1) and give an explanation of the origin of the obs erve d BAU 



(Sec. 4.2 1 consistently with the G constraint and the low energy neutrino data (Sec. 



4.3 



4.1 The Inflaton's Decay 

When nMHI is over, the inflaton continues to roll down towards the SUSY vacuum, Eq. (9). 
There is a brief stage of tachyonic preheating d40t which does not lead to significant particle 
production fill . Soon after, the inflaton settles into a phase of damped oscillations initially 
around zero - where Vhio has a maximum - and then around one of the minima of Vhio- 
Whenever the inflaton passes through zero, particle production may occur creating mostly 
superheavy bosons via the mechanism of instant preheating ( 42 1. This process becomes more 
efficient as A decreases, and further numerical investigation is required in order to check the 
viability of the non-thermal leptogenesis scenario for small values of A. For this reason, we 
restrict to A's larger than 0.001, which ensures a less frequent passage of the inflaton through 
zero, weakening thereby the effects from instant preheating and other parametric resonance 
effects - see Appendix B of Ref. (5 1. Intuitively the reason is that larger A's require larger c^'s, 
see Eq. (31) , diminishing therefore hf given by Eq. (29) , which sets the amplitude of the very 
first oscillations. 

Nonetheless the standard perturbative approach to the inflaton decay provides a very efficient 
decay rate. Namely, at the SUSY vacuum v c H and v c H acquire the v.e.vs shown in Eq. (9) 
giving rise to the masses of the (canonically normalized) inflaton Sh = (h — 2Mpg) / /o and 
RH neutrinos, xK, which are given, respectively, by 

^AMpc Xi v cMk Q 
(a) m Y = Vijjrf and (b) M iv , = 2 ' v ^ (33) 
Wfo MsV/o 



where /o is defined below Eq. 25b and fa = fo + 24c^?-rcL. ~ /q. Here, we assume the 
existence of a term similar to the i 



■ second one inside In of Eq. |7al for v? too. 

For larger A's (/) = }(h = 2Mps) ranges from 3 to 90 and so m\ is kept independent of A and 
almost constant at the level of 10 13 GeV. Indeed, if we express Sh as a function of Sh through 
the relation 

^ ~ J where J = J 1 + \m\f\ ((h)) = ^\+2\c\m\ s (34) 
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we find 

^ v^Mps ^ J*r_ ^ lO^mp ^ 3 1Ql3 ^ for A 1Q-4 ^ u io _ 3 
/o/o 2V3c K 4.2^3 4.2V 6mpg 

where we make use of Eq. j3T) - note that /q ~ 1. The derivation of the (s)particle spectrum, 
listed in Table 2, at the SUSY vaccum of the model reveals (5i that perturbative decays of 
Sh into these massive particles are kinematically forbidden and therefore, narrow parametric 
resonance (40 1 effects are absent. Also Sh can not decay via renormalizable interaction terms 
to SM particles. 

The inflaton can decay into a pair of vf's through the following lagrangian terms: 

Clv < = " A,Vl mT f " 12c7 ^ s ) + hx - ' (36) 

From Eq. |36| we deduce that the decay of Sh into vf is induced by two lagrangian terms. The 
first one originates exclusively from the non-renormalizable term of Eq. (5) - as in the case of a 
similar model in Ref. (|30|. The second term is a higher order decay channel due to the SUGRA 
lagrangian - cf. Ref. (43 1. The interaction in Eq. (36) gives rise to the following decay width 



c 2 



4M 2 ~ M~ f 3/2 

1 1 = 64^ 1 " With ^ = m£ V ( 2 " 12C7 ^s) • ^ 



where Mp: is the Majorana mass of the vj's into which the inflaton can decay. The 



implementation - see Sec. 4.3 - of the seesaw mechanism for the derivation of the 



light-neutrinos masses, in conjunction with the Gpg prediction m 3D ~ nit and our assumption 



that < f«2D <C ff*3D - see Sec. 5.1 - results to 2M 3 j7i > mj. Therefore, the kinematically 
allowed decay channels of Sh are those into v°- with j = 1 and 2. Note that the decay of the 

inflaton to the heaviest of the Vj's (V3) is also disfavored by the G constraint - see below. 

In addition, there are SUGRA-induced (43 1 - i.e., even without direct superpotential couplings 
- decay channels of the inflaton to the MSSM particles via non-renormalizable interaction 
terms. For a typical trilinear superpotential term of the form W y = yXYZ, we obtain the 
effective interactions described by the langrangian part 

r3 /2 

Ay = 6yc TC rc/o f h fx{p Y {p z + 9{j, x {p z + Zf x f Y ) + hx. , (38) 

Mp 2/o V / 

where y is a Yukawa coupling constant and ip x , tpy an d tpz ar e the chiral fermions associated 
with the superfields X, Y and Z. Their scalar components are denoted with the superfield 
symbol. Taking into account the terms of Eq. |2) and the fact that the adopted SUSY GUT 
predicts YU for the 3rd generation at Mpg, we conclude that the interaction above gives rise 
to the following 3-body decay width 

r * - = % - <39) 
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with 1/33 — (0.55 — 0.7) being the common Yukawa coupling constant of the third generation 
computed at the m\ scale, and summation is taken over color, weak and hypercharge degrees 
of freedom, in conjunction with the assumption that m\ < IMygc. 

Since the decay width of the produced is much larger than Tj- see below - the reheating 
temperature, T^, is exclusively determined by the inflaton decay and is given by (1441 



«rh 



72 xl/4 



5n 2 g, 



V T i m P with ^ = T mc + T m c + Tiy, (40) 



where counts the effective number of relativistic degrees of freedom at temperature T^. 
For the MSSM spectrum plus the particle content of the superfields P and P we find g* ~ 
228.75 + 4(1 + 7/8) = 236.25. 

4.2 Lepton-Number and Gravitino Abundances 

If T r i, <C Mffc, the out-of -equilibrium condition (2) for the implementation of nTL is 
automatically satisfied. Subsequently decay into H u and L* via the tree-level couplings 
derived from the second term in the RHS of Eq. J5J. Interference between tree-level and 
one-loop diagrams generates a lepton-number asymmetry (per Vj decay) £j {2), when CP 
conservation is violated. The resulting lepton-number asymmetry after reheating can be 
partially converted through sphaleron effects into baryon-number asymmetry. In particular, 
the B yield can be computed as 

(a) Y B = -0.35Y L with (b) Y L =l\ ^ £ (41) 



The numerical factor in the RHS of Eq. j4Tp ) comes from the sphaleron effects, whereas the one 
(5/4) in the RHS of Eq. (ilfr ) is due to the slightly different calculation {2} of T rh - cf. Ref. ©. 



In the major part of our allowed parameter space - see Sec. 5.2 - Y\ ~ Y\y and so the involved 
branching ratio of the produced is given by 

rW + Tgg ^ Tgg _ n 2 (l- 12c n m 2 FS ) 2 M 2 2 ~ c 
Ti '" T Iy 72c|y2 3m 4 g m 2 

For M 2 vf — (10 11 — 10 12 ) GeV the ratio above takes adequately large values so that Yj_ is 
sizable. Therefore, the presence of more than one inflaton decay channels does not invalidate 
the scenario of nTL. 

It is worth emphasizing, however, that if M\ v c < 10T r j„ part of the Yi can be washed out 
due to v\ mediated inverse decays and AL = 1 scatterings - this possibility is analyzed in 
Ref. (27i. Trying to avoid the relevant computational complications we limit ourselves to cases 
with M\-yc > lOTj-h, so as any washout of the non-thermally produced Yi is evaded. On the 
other hand, Yi is not erased by the AL = 2 scattering processes 1 45 1 at all temperatures T with 
100 GeV < T < T lh since Y L is automatically protected by SUSY {46) for 10 7 GeV < T < T rh 
and for T < 10 7 GeV these processes are well out of equilibrium provided that that mass of 
the heaviest light neutrino is 10 eV. This constraint, however, is overshadowed by a more 



stringent one induced by WMAP7 data 1 16 1 - see Sec. 5.1 
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The required for successful nTL T r ^ must be compatible with constraints on the G abundance, 
Yg, at the onset of nucleosynthesis (BBN). This is estimated to be 1 19): 

Yg ~ cgT rh with eg = 1.9 ■ lCT 22 /GeV, (43) 

where we assume that G is much heavier than the gauginos. Let us note that non-thermal 
G production within SUGRA is (43) also possible. However, we here prefer to adopt the 
conservative approach based on the estimation of Yg via Eq. |43| since the latter G production 
depends on the mechanism of SUSY breaking. 



Both Eqs. 1 41 and (43) yield the correct values of the B and G abundances provided that no 
entropy production occurs for T < T rh . This fact can be easily achieved within our setting. 
The mass spectrum of the P-P system is comprised by axion and saxion P_ = (P — P) / yfl, 
axino = (fp — ip P )/ \fl, a higgs, P+ = (P + P) / y/2, and a higgsino, xp + = (ipp +ipp)/ V2, 
with mass of order 1 TeV and tp denoting a Weyl spinor. The higgs and higgsinos can decay 
to lighter higgs and higgsinos before domination (36). Regarding the saxion, P_, we can 
assume that its decay mode to axions is suppressed (w.r.t the ones to gluons, higgses and 
higgsinos 1 47 : 48 1) and the initial amplitude of its oscillations is equal to /« ~ 10 12 GeV. Under 
these circumstances, it can (47) decay before domination too, and evades (48 1 the constraints 
from the effective number of neutrinos for the f a 's and T rh 's encountered in our model. As a 
consequence of its relatively large decay temperature, the LSPs produced by the saxion decay 
are likely to be thermalized and therefore, no upper bound on the saxion abundance is (48) to 
be imposed. Finally, axino can not play the role of LSP due to its large expected mass and the 
relatively high Trf/s encountered in our set-up which result to a large Cold Dark Matter (CDM) 
abundance. Nonetheless, it may enhance non-thermally the abundance of a higgsino-like 
neutralino-LSP, rendering it a successful CDM candidate. 

4.3 Lepton-Number Asymmetry and Neutrino Masses 

As mentioned above, the decay of v\ and emerging from the Sh decay, can generate a lepton 
asymmetry, e, (with i = 1,2) caused by the interference between the tree and one-loop decay 
diagrams, provided that a CP-violation occurs in h N {j's. The produced e, can be expressed in 
terms of the Dirac mass matrix of v,-, mrj, defined in a basis (called vf-basis henceforth) where 
V c . are mass eigenstates, as follows: 

Im [(m^rno) 2 . / \ 
^ 87r(H„) 2 (mJ ) m D ) !i - V v 'J 

where we take (H u ) ~ 174 GeV, for large tan & and 

M ]v< , T, V r {m f D m D )i, 

X " = md * : = = 8^F (44b) 

(with i,j = 1,2,3). Also Fy an d Fs represent, respectively, the contributions from vertex and 
self-energy diagrams which in SUSY theories read (49) 



Fy (x) = -iln (l + x- 2 ) and F s (x,y,z) = — (44c) 
v ' (x 2 - 1) + (x 2 z - y) 
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Note that for strongly hierarchical M^ 's with x» ~> 1 and x,y 3> J/;,1/;, we obtain the 
well-known approximate result i26ll27l 

fv + F S --3/4. (45) 



The involved in Eq. j44af can be diagonalized if we define a basis - called weak basis 
henceforth - in which the lepton Yukawa couplings and the SU(2)l interactions are diagonal 
in the space of generations. In particular we have 



(J + m D !i c+ = d D = diag(m 1D ,M 2D ,m 3D ) 



(46) 



where U and U c are 3x3 unitary matrices which relate L, and vj (in the v ( c -basis) with the 
ones L; and vf in the weak basis as follows: 



l! = LU and v cl = U c v c 



(47) 



Here, we write LH lepton superfields, i.e. SU (2)l doublet leptons, as row 3-vectors in family 
space and RH anti-lepton superfields, i.e. S(I(2) L sing let a nti-leptons, as column 3-vectors. 
Consequently, the combination ni^niT) appeared in Eq. 1 44a 1 turns out to be a function just of 
drj and U c . Namely, 

m^m u = U^d^dull . (48) 

The connection of the nTL scenario with the low energy neutrino data can be achieved through 
the seesaw formula, which gives the light-neutrino mass matrix m v in terms of niffj and Mfpc. 
Working in the v ( c -basis, we have 



ty = —niD d v } wIq, where d v c = diag (M lv - C , M 2 pr,M ; 



3f" 



(49) 



with Mjpr < M 2 yc < M 3l 7c real and positive. Solving Eq. |46j w.r.t m D and inserting the 
resulting expression in Eq. |49| we extract the mass matrix 



U f mM* 



-dvU'd-^U^dv, 



which can be diagonalized by the unitary PMNS matrix satisfying 

m v = U* diag (?«!„, M 2l/ ,ff%,) ilj 

and parameterized as follows: 



Uy = 



with c,. 



C12C13 



s 12 c 13 



Sl3 e 



-C23S12 - S 2 3C\2 s 13 e ' S c 23 c 12 ~ S23Sl2Sl3 e!<5 S 23 C 13 



, S23S12 - C23 c 12Sl3£ 



-s 23 c 12 - c 23 s 12 s 13 e c 23 c 13 _ 



-ifi/2 
e-*>2/2 
1 



(50) 



(51) 



(52) 



cos 8 ij, Sij := sindij, S the CP-violating Dirac phase and q>\ and <p 2 the two 
CP-violating Majorana phases. 

Following a bottom-up approach, along the lines of Ref. i26H28l , we can find m v via Eq. |5T| 
using as input parameters the low energy neutrino observables, the CP violating phases and 
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adopting the normal or inverted hierarchical scheme of neutrino masses. Taking also m,rj as 
input parameters we can construct the complex symmetric matrix 

W = -d^rfivd^ 1 = U c d v cU cJ (53) 

- see Eq. (50) - from which we can extract d v c as follows: 

d~ 2 = U ci WW f U c . (54) 

Note that WW f is a 3 x 3 complex, hermitian matrix and can be diagonalized following the 
algorithm described in Ref. l50l . Having determined the elements of U c and the Mjp/s we can 
compute m D through Eq. |48j and the £,'s through Eq. l |44a} . 



5. Constraining the Model Parameters 



We exhibit the constraints that we impose on our cosmological set-up in Sec. 5.1 and delineate 
the allowed parameter space of our model in Sec.|5.2| 



5.1 Imposed Constraints 

The parameters of our model can be restricted once we impose the following requirements: 

1 . According to the inflationary paradigm, the horizon and flatness problems of the standard 
Big Bang cosmology can be successfully resolved provided that the number of e-foldings, 
N„, that the scale fc* = 0.002/Mpc suffers during nMHI takes a certain value, which 
depends on the details of the cosmological model. Employing standard methods JBTl , we 
can easily derive the required N* for our model, consistently with the fact that the P — P 
system remains subdominant during the post-inflationary era. Namely we obtain 

2. The inflationary observables derived in Sec. |3.2| are to be consistent with the fitting Jl6t of 
the WMAP7, BAO and Hq data. As usual, we adopt the central value of An, whereas we 
allow the remaining quantities to vary within the 95% confidence level (c.l.) ranges. Namely, 

(a) h-ji ~ 4.93 ■ 1CT 5 , (b) n s = 0.968 ± 0.024, (c) - 0.062 < a s < 0.018 and (d) r < 0.24 

(56) 

3. The scale Mpg can be determined by requiring that the v.e.vs of the Higgs fields take the 
values dictated by the unification of the gauge couplings within the MSSM. As we now 
recognize - cf. Ref. (5 1 - the unification scale Mqjt — 2 ■ 10 16 GeV is to be identified with 
the lowest mass scale of the model in the SUSY vacuum, Eq. |9), in order to avoid any extra 
contribution to the running of the MSSM gauge couplings, i.e., 

(a) 8 ~7W = MGUT = ^ (b) ^ = (57) 

The requirement 2c™ x — c^j > sets an upper bound c-ji < 2c™ ax ~ 1.8 ■ 10 3 , which 



however can be significantly lowered if we combine Eqs. (55) and 128) - see Sec. 5.2.1 
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4. For the realization of nMHI , we assume that takes relatively large values - see 
e.g. Eq. (17) . This assumption may ((52) [53} jeopardize the validity of the classical 
approximation, on which the analysis of the inflationary behavior is based. To avoid 
this inconsistency - which is rather questionable JTTJ |54) though - we have to check 
the hierarchy between the ultraviolet cut-off, A = nip/c-n, of the effective theory and 
the inflationary scale, which is represented by Vhi^*) 1 ^ 4 or, less restrictively by the 
corresponding Hubble parameter, fi* = Vfji(^*) 1 / \/3mp. In particular, the validity of 
the effective theory implies i52ll53l 

(a) V H i(M 1/4 < A or (b) H* < A for (c) c n > 1. (58) 

5. As discussed in Sec. |4.2| to avoid any erasure of the produced Yi and to ensure that the 
inflaton decay to v 2 is kinematically allowed we have to bound Mjpc and M 2 j?r respectively 
as follows: 

(a) A%* > 10T rh and (b) m Y > 2M 2 p => M 2 yc < A "! F ~ 1.5 ■ 10 13 GeV, (59) 

where we make use of Eq. |35) . Recall that we impose also the restriction A > 0.001 which 
allows us to ignore effects of instant preheating (5"ll42l. 

6. As discussed below Eq. |2|, the adopted GUT predicts YU at Mqjjj. Assuming negligible 
running of m 3D from Mquj until the scale of nTL, Aj_, which is taken to be A^ = m t , we 
end up with the requirement: 

m 3D (m I ) = m,(mi) ~ (100 - 120) GeV. (60) 

where m t is the top quark mass and the numerical values correspond to J/33(»2i) = (0.55 — 

0. 7) - cf. Ref. (55) - found (32: 561 working in the context of several MSSM versions with 
tan /3 ~ 50 and taking into account the SUSY threshold corrections. As regards the lighter 
generation, we limit ourselves in imposing just a mild hierarchy between mjrj and m 2 D, 

1. e., niio < W2D m 3D since it is not possible to achieve a simultaneous fulfilment of all 
the residual constraints if we impose relations similar to Eq. {60} - cf. Ref. f2"5H2"7t . 

7. From the solar, atmospheric, accelerator and reactor neutrino experiments we take into 
account the inputs listed in Table [3] on the neutrino mass-squared differences Am^ and 
Am^, on the mixing angles and on the CP-violating Dirac phase, S for normal [inverted] 
neutrino mass hierarchy (23} - see also Ref. (24 1. In particular, m, v 's can be determined via 
the relations: 

\l m \v + ^ m 3i' for normally ordered (NO) m v 's 

or (61) 
y m\ v + I Am^j |, for invertedly ordered (IO) m v 's 

The sum of m iv 's can be bounded from above by the WMAP7 data dl6t 

Eim iv < 0.58 eV (62) 

at 95% c.l. This is more restrictive than the 95% c.l. upper bound arising from the effective 
electron neutrino mass in /J-decay f57): 



m 2v 



m\ v + Am|j and 



m 3 v 



miv 



(63) 
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Parameter 


Best Fit ±la 


Normal 


Inverted 


Hierarchy 


Am^/io-SeV 2 
Aw| x /10- 3 eV 2 


7.62 ±0.19 


9 co+0.08 

-o.io 


9 yi+0.10 
z - -0.07 


sin 2 9i2 
sin 2 13 
sin 2 023 


320+ 0015 


026+ - 003 
0.49t°;°| 


n 027 +0.003 

-0.004 

53+ 05 


S/n 


0JB$* 


0.07 



Table 3. Low energy experimental neutrino data for normal or inverted hierarchical neutrino masses. In 
the second case the full range (0 — In) is allowed at lu for the phase S. 

However, in the future, the KATRIN experiment (581 expects to reach the sensitivity of 
mp ~ 0.2 eV at 90% c.l. 

8. The interpretation of BAU through nTL dictates QU at 95% c.l. 



Y B = (8.74 ± 0.42) ■ 10 



-11 



32 < 10 n Y B < 9.16. 



(64) 



9. In order to avoid spoiling the success of the BBN, an upper bound on Yg is to be imposed 
depending on the G mass, Mg, and the dominant G decay mode. For the conservative case 
where G decays with a tiny hadronic branching ratio, we have H9\ 



10 -14 

for mg 

io- 12 



0.69 TeV 
10.6 TeV 
13.5 TeV. 



(65) 



As we see below, this bound is achievable within our model only for mg > 10 TeV. Taking 
into account that the soft masses of the scalars are not necessarily equal to nig, we do not 
consider such a restriction as a very severe tuning of the SUSY parameter space. Using 
Eq. |43| the bounds on Yg can be translated into bounds on T r j v Specifically we take 
T rh ~ (0.53 - 5.3) ■ 10 s GeV [T rh ~ (0.53 - 5.3) ■ 10 9 GeV] for Yg ~ (0.1 - 1) ■ 10~ 13 
[Yg ~ (0.1 - 1) ■ 10" 12 ]. 



Let us, finally, comment on the axion isocurvature perturbations generated in our model. 
Indeed, since the PQ symmetry is broken during nMHI, the axion acquires quantum 
fluctuations as all the almost massless degrees of freedom. At the QCD phase transition, 
these fluctuations turn into isocurvature perturbations in the axion energy density, which 
means that the partial curvature perturbation in axions is different than the one in photons. 
The results of WMAP put stringent bounds on the possible CDM isocurvature perturbation. 



5 Constraining the Model Parameters 



19 



Namely, taking into account the WMAP7, BAO and Hq data on the parameter kq we find the 
following bound for the amplitude of the CDM isocurvature perturbation 



\S C \ = h nx l-^— < 1.5 ■ 1CT 5 at 95% c.l. (66) 



On the other, |5 C | due to axion, can be estimated by 

| 5c l = 9± Jk^ with £±~a( f' — V' 175 (67) 

where O fl [Q c ] is the axion [CDM] density parameter, (pp* ~ 10 16 GeV d36t denotes the field 
value of the PQ scalar when the cosmological scales exit the horizon and 0\ is the initial 



misalignment angle which lies (36) in the interval [—n/6,n/6\. Satisfying Eq. '66' requires 
\&l\ < n/70 which is a rather low but not unacceptable value. Therefore, a large axion 
contribution to CDM is disfavored within our model. 

5.2 Numerical Results 

As can be seen from the relevant expressions in Secs.|5]and|4] our cosmological set-up depends 
on the parameters: 

A, Ah, A h , k s , g, i/33, m lvi m iD , q>i and q> 2 , 

where m( v is the low scale mass of the lightest of i/;'s and can be identified with m\ v [m^ v \ for 
NO [IO] neutrino mass spectrum. Recall that we determine Mpg via Eq. (57) with g = 0.7. We 
do not consider c-ji and A, v c as independent parameters since c-ji is related to m via Eq. (3T) 
while Ki v c can be derived from the last six parameters above which affect exclusively the Y L 
calculation and can be constrained through the requirements 5 - 9 of Sec. |5.1| Note that the 
A, v c's can be replaced by M^/s given in Eq. \33\ >) keeping in mind that perturbativity requires 
Ajyc < V / 47r or M- t yc < 10 16 GeV. Note that if we replace Mg with my in Eq. ji}, we obtain a 
tighter bound, i.e., M;p= < 2.3 • 10 15 GeV. Our results are essentially independent of Ah, A^ 
and k s , provided that we choose some relatively large values for these so as m~_, m~ and m| 
in Tableware positive for A < 1. We therefore set A H = A H = 0.5 and kg = 1 throughout our 
calculation. Finally T r ^ can be calculated self-consistently in our model as a function of m\, 



M-2yc S> Mjft and the unified Yukawa coupling constant 1/33 - see Sec. 4.1 - for which we take 
y 33 = 0.6. 

Summarizing, we set throughout our calculation: 

k s = 1, Ah = Ah = 0.5, g = 0.7 and y 33 = 0.6. (68) 

The selected values for the above quantities give us a wide and natural allowed region for the 
remaining fundamental parameters of our model, as we show below concentrating separately 
in the inflationary period (Sec. 5.2.1} and in the stage of nTL (Sec. 5.2.2} . 



5.2.1 The Inflationary Stage 

In this part of our numerical code, we use as input parameters ?/*,m 2 D Wip and en- For 
every chosen c-ji > \ and tn 2 u, we restrict A and h* so that the conditions Eq. (55) and (56^) 
are satisfied. In our numerical calculations, we use the complete formulas for the slow-roll 
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10 




T t /10°GeV 



10 1 - 



10 " 



(a) 



10" 




Figure 1 . The allowed (by all the imposed constraints) values ofc^ (solid line), - given by Eq. \4(fy 
(dashed line), m\ (dot-dashed line) and M 2 vc (dotted line) [hf (solid line) and h t (dashed line)] versus A 
(a) [(b)] forks = 1, Ah = A H = 0.5 and 1/33 = 0.6. The light gray and gray segments denote values of the 
various quantities satisfying Eq. j5Sb) too, whereas along the light gray segments we obtain h t > rap. 



parameters and An in Eqs. |25a| , j25b} and (30) and not the approximate relations listed in 
Sec. |3.2| for the sake of presentation. Our results are displayed in Fig. [lj where we draw the 
allowed values of c-ji (solid line), T r h (dashed line), the inflaton mass, mj (dot-dashed line) and 
M-2vc (dotted line) - see S ec.|4.1| - [hf (solid line) and h* (dashed line)] versus A (a) [(b)] for the 
m 2 D's required from Eq. 
Fig. [2]- see Sec. 5 

The lower bound of the depicted lines comes from the saturation of the Eq. |58t). The 



|64| and for the parameters adopted along the black dashed line of 
The required via Eq. 1 55 remains almost constant and close to 54.5. 



constraint of Eq. (58p ) is satisfied along the various curves whereas Eq. j58p ) is valid only 
along the gray and light gray segments of these. Along the light gray segments, though, we 
obtain /i* > nip. The latter regions of parameter space are not necessarily excluded, since 
the energy density of the inflaton remains sub-Planckian and so, corrections from quantum 
gravity can still be assumed to be small. As c-ji increases beyond 906, fo becomes much larger 
than 1, derived by Eq. |28| starts decreasing and therefore, nMHI fails to fulfil Eq. 1 55 1. 
This can be understood by the observation that N„, approximated fairly by Eq. 1 29 >, becomes 
monotonically decreasing function of c-ji for c-ji > c^ ax where c™" can be found by the 
condition 

dN^^ 3hl {c^ x -cn 



4m| 



-max 







(69) 



where c^ ax is defined in Eq. (57p) and Eq. 1 57 i) is also taken into account. As a consequence, 



SY GUT provides us with a clear upper bound of c-ji. All in 

(70) 



the embedding of nMHI in a 
all, we obtain 

0-001 < A < 0.042 and 1 < c n < 907 for N* ~ 54.5 
When c-ji ranges within its allowed region, we take M PS ~ (2.87 - 4) ■ 10 16 GeV. 

From Fig. [l]-(a), we can verify our analytical estimation in Eq. pi") according to which A is 
proportional to c-ji- On the other hand, the variation of hf and h* as a function of c-ji - drawn 
in Fig.[l]-(b) - is consistent with Eqs. {27} and (29| . Letting A or c-ji vary within its allowed 
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region in Eq. (70) , we obtain 

n s ~ 0.964, a s ~ -6.35 ■ 10~ 4 and r ~ 3.6 ■ 10~ 3 . (71) 

Clearly, the predicted a s and r lie within the allowed ranges given in Eq. (56j)) and Eq. (56J:) 
respectively, whereas n s turns out to be impressively close to its central observationally 
favored value - see Eq. (56fc ) and cf. Ref. l lT2l . 

From Fig.IIWa) we can conclude that mi is kept independent of A and almost constant at the 
level of 10" GeV, as anticipated in Eq. (35). From the same plot we also remark that for 



A < 0.03, T r h remains almost constant since I\y dominates over r I2 j;c and /q ~ 1 - see Eq. 1 39 
For A > 0.03, /q ~ 1 + Hc-jim^ starts to deviate from unity and so, T r j, increases with c-ji or 
A as shown in Fig. [I] The required by Eq. (64) Mq$c follows the behavior of the required m^o 
- see Fig.[2)(a) of Sec. [5^2| 

5.2.2 The Stage of non-Thermal Leptogenesis 

In this part of our numerical program, for a given neutrino mass scheme, we take as input 
parameters: mgy, m, D , q>\, f2 and the best-fit values of the neutrino parameters listed in Table|3] 
We then find the renormalization group (RG) evolved values of these parameters at the scale of 
nTL, Ai, which is taken to be Aj_ = m.\, integrating numerically the complete expressions of 
the RG equations - given in Ref. (29) - for m,„, 0«, S, q>\ and q>2- In doing this, we consider the 
MSSM with tan/3 ~ 50, favored by the preliminary LHC results - see, e.g., Ref. 132; 56 1 - as 
an effective theory between and a SUSY-breaking scale, MgusY — 1-5 TeV. Following the 
procedure described in Sec. |4.3| we evaluate M,^ at Aj,. We do not consider the running of 
m (D and Mj$c and therefore, we give their values at A L . 

We start the exposition of our results arranging in Table [4] some representative values of the 
parameters leading to the correct BAU for A = 0.01 and c-ji = 220 and normally hierarchical 
(cases A and B), degenerate (cases C, D and E) and invertedly hierarchical (cases F and G) 
m v 's. For comparison we display the B-yield with (Yg) or without (Yj^) taking into account 
the RG effects. We observe that the two results are more or less close with each other. In all 
cases the current limit of Eq. (62) is safely met - the case D approaches it -, while mp turns out 
to be well below the projecteasensitivity of KATRIN ( 58 1. Shown are also the obtained T rh 's, 
which are close to 10 9 GeV in all cases, and the corresponding Yg's, which are consistent with 



Eq. 1 65 



for »xg > 11 TeV. 



From Table[4]we also remark that the achievement of Y B within the range of Eq. (64) dictates a 
clear hierarchy between the M-ffc 's, which follows the imposed hierarchy in the sector of m, D 's 
- see paragraph 6 of Sec. |5.1| This is expected since, in the limit of hierarchical m;p's, the M,-pc's 
can be approximated by the following expressions I 




for NO m v 's 

(M lpc ,M 2 ^A%0~<( ( V X 12 l^ "rL\ . TO , (72a) 

for IO m v s 

Indeed, we see e.g. that for fixed /, the Mac's depends exclusively on the trip's and Mjgc 
increases when ni( v decreases with fixed m 3D . As a consequence, satisfying Eq. (59fc ) pushes 
the nijrj's well above the mass of the quark of the first generation. Similarly, the /«2d' s required 
by Eq. (64) turns out to be heavier than the quark of the second generation. Also, the required 
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Parameters 


Cases 




A 


B 


C 


D 


E 


F 


G 




Normal 


Degenerate 


Inverted 




Hierarchy 


Masses 


Hierarchy 


Low Scale Parameters 


m lv /0.1 eV 


0.01 


0.1 


0.5 


1. 


0.7 


0.5 


0.49 


m 2v /0.1 eV 


0.088 


0.13 


0.5 


1. 


0.7 


0.51 


0.5 


m 3v /0.1 eV 


0.5 


0.5 


0.71 


1.1 


0.5 


0.1 


0.05 


E;m,v/0.1eV 


0.6 


0.74 


1.7 


3.1 


1.9 


1.1 


1 


mp/0.1 eV 


0.03 


0.013 


0.14 


0.68 


0.45 


0.33 


0.4 


n 


n/2 


n/2 


n/2 


-n/2 


-n/2 


-n/2 


-7T/4 


(p2 





-n/2 


-n/2 


n 


n 


n 





Leptogenesis-Scale Parameters 


m 1D /GeV 


0.4 


0.3 


0.8 


1 


0.9 


0.9 


0.9 


m 2 D/GeV 


9.2 


3 


6.5 


8.6 


3.95 


6.6 


9.2 


m 3D /GeV 


120 


100 


100 


120 


120 


110 


110 


M li;c /10 10 GeV 


4.4 


4.7 


3.6 


1.2 


1.5 


1.6 


1.7 


M 2V c /10 12 GeV 


2.7 


0.6 


0.65 


1.5 


0.9 


1.5 


2.8 


M 357 c/10 14 GeV 


27 


0.28 


0.46 


0.4 


0.4 


3.8 


10 


Resulting B-Yield 


10 n Y° 


8.75 


8.9 


8.6 


8.63 


8.9 


9. 


8.76 


10 U Y B 


9.3 


8.7 


8.5 


8.4 


9.7 


8.8 


9.2 


Resulting T rh and Y^ 


T rh /10 9 GeV 


1.2 


0.89 


0.89 


0.99 


0.91 


0.99 


1.2 


lO^Yg 


2.4 


1.7 


1.7 


1.9 


1.7 


1.88 


1.88 



Table 4. Parameters yielding the correct BAU for various neutrino mass schemes for A = 0.01 and 
c K = 220. 

by seesaw M^yc 's are lower in the case of degenerate Vj spectra and can be as low as 3 • 10 13 GeV 
in sharp contrast to our findings in Ref . 1 5 1, where much larger M^yc 's are necessitated. An 
order of magnitude estimation for the derived Ei's can be achieved by 1 25 : 26 1 

3M 2l7l f(^) forNOm v 's 
on{M u ) ^ mv fcj-iOwVs 

which is rather accurate, especially in the case of IO m v 's. 
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q 8 



i 1 


1 • 

i \ : 


m 1v / eV 


m 1D / GeV 


m 3D /GeV 


% 


% 


I 1 


i \ 


0.001 


0.4 


120 


nil 







i ■, 

V \ " 


0.01 


0.3 


100 


nil 


1,12 


/ 1 

1 


\ '\ 


0.05 


0.8 


100 


n!2 


-n/2 


/ 1 


\ 


0.1 


1 


120 


-n/2 


71 


/ 1 


\ 












,•' / 
/ 


■ 


m 3v /eV 


m 1D / GeV 


m 3D /GeV 


<P, 


% 




0.05 


0.9 


120 


-u/2 


n 






0.01 


0.9 


110 


-it/2 


n 






0.005 


0.9 


110 


-it/4 






10" 



m lv (eV) 


A (10- 2 ) 


m 2D (GeV) 


Mu;c(10 lu GeV) 


M 2V c{\Q 12 GeV) 


A% c (10 14 GeV) 


0.001 


0.1 - 1.8 


5.5 - 17 


4.7 


1-9.8 


32 




2.4 - 4.1 


17-10 


4.7 


17-10 


32 


0.01 


0.1-2 


1.5-13 


3.5-6 


0.26 - 8.75 


0.33 - 0.4 




2-4.1 


15 - 3.4 


6-5.3 


11 - 0.84 


0.42 - 0.34 


0.05 


0.1-2. 


3-28 


3.5-4 


2.1-11 


0.5 - 0.6 




2.1 - 4.1 


28 - 7.3 


4-3.9 


11-0.8 


0.6 - 0.5 


0.1 


0.1-1.9 


4-23 


1.3 


0.4-8 


0.4 - 0.5 




2.2-4.1 


23-9 


1.3 


10-2 


0.5 - 0.4 


f«3v(eV) 


A (W- z ) 


m 2D (GeV) 


M lt7 c(10 iu GeV) 


M 2l7c (10 12 GeV) 


M 3i7c (10 14 GeV) 


0.05 


0.1-2 


2.1 - 12 


1.5-1.6 


0.28 - 7.8 


0.48 - 0.56 




2.1 - 4.1 


13 - 4.4 


1.6 


8.7-1.2 


0.57 - 0.49 


0.01 


0.1-2 


2.8 - 17 


2.2 


0.4-11 


5.4-5.6 




2-4.1 


17-7.3 


2.2 


11-2 


5.6-5.4 


0.005 


0.1 - 1.8 


5-17 


1.8 


0.7-10.1 


12 




1.8-4.1 


17.7 - 10 


1.8 


10.5 - 3.5 


12 



Figure 2. Contours on the A — m^o plane, yielding the central Yg in Eq. \64\ , consistently with the 
inflationary requirements, for Ah = Afj = 0.5,^5 = 1 andt/33 = 0.6 and various (mi v ,mYQ,q>\,q>7)'s 



indicated next to the graph (c) and NO [IO] m„ 
are also shown in the table included. 



's (black [gray] lines). The corresponding ranges of M^c 's 
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To highlight further our conclusions inferred from Table[4] we can fix m^ v (m lv for NO m,„'s 
or m^y for IO mj v 's) ffxjrj, m^, <P\ and <P2 to their values shown in this table and vary m^o 
so that the central value of Eq. (64) is achieved. This is doable since, according Eq. 1 172a} , 
variation of m 2 -Q induces an exclusive variation to M 2 yc which, in turn, heavily influences ej, — 
see Eqs. (44a) and (45) - and Yi - see Eqs. (IT) and (42). T he resulting contours in the A — m 2 D 
plane are presented in Fig. [5]- since the range of Eq7|64) is very narrow the possible variation 
of the drawn lines is negligible. The resulting M^/s are displayed in the table included. The 
conventions adopted for the types and the color of the various lines are also described next 
to the graph (c) of Fig. [3] In particular, we use black [gray] lines for NO [IO] m,v's. The black 
dashed and the solid gray line terminate at the values of m 2 u beyond which Eq. (64) is non 
fulfilled due to the violation of Eq. (59b). 



In all cases, two disconnected allowed domains arise according to which of the two 



contributions in Eq. 1 36 dominates. The critical point (A c , cn c ) is extracted from: 



1 - \2c Uc ml s = 0^c nc = c^ ax /2 ~ 453 or A c ~ 10" 4 c^ ax /4.2 ~ 0.021 (73) 

where we make use of Eq . (57) and Eq. (31) in the intermediate and the last step respectively. 
From Eqs. (40) , (41) and (42) one can deduce that for A < A c , T rh remains almost constant; 
T^/Tj decreases as c-ji increases and so the M 2 yc's, which satisfy Eq. (64) , increase. On 
the contrary, for A > A c , T\ 2 yc /Tj is independent of c-ji but T r ^ increases with c-ji and so the 
fulfilling Eq. (64) M 2 y C 's decrease. 

Summarizing, we conclude that our scenario prefers the following ranges for the Mg/c's: 

1 < /10 10 GeV < 6, 0.6 < M 2t7l /10 12 GeV < 20, 0.3 < M 3pi /10 14 GeV < 30, (74a) 
while the m\x> ar *d wi2D are restricted in the ranges: 

0-3 < m 1D /GeV < 1, 1.5 < m 2D /GeV < 20. (74b) 



6. Conclusions 

We investigated the implementation of nTL within a realistic GUT, based on the PS gauge 
group. Leptogenesis follows a stage of nMHI driven by the radial component of the Higgs 
field, which leads to the spontaneous breaking of the PS gauge group to the SM one with 
the GUT breaking v.e.v identified with the SUSY GUT scale and without overproduction of 
monopoles. The model possesses also a resolution to the strong CP and the ]i problems of the 
MSSM via a PQ symmetry which is broken during nMHI and afterwards. As a consequence 
the axion cannot be the dominant component of CDM, due to the present bounds on the axion 
isocurvature fluctuation. Moreover, we briefly discussed scenaria in which the potential axino 
and saxion overproduction problems can be avoided. 

Inflation is followed by a reheating phase, during which the inflaton can decay into the 
lightest, v\, and the next-to-lightest, v\, RH neutrinos allowing, thereby for nTL to occur 
via the subsequent decay of v\ and vS. Although other decay channels to the SM particles 
via non-renormalizable interactions are also activated, we showed that the production of the 
required by the observations BAU can be reconciled with the observational constraints on 
the inflationary observables and the G abundance, provided that the (unstable) G masses are 
greater than 11 TeV. The required by the observations BAU can become consistent with the 
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present low energy neutrino data, the restriction on m 3D due to the PS gauge group and the 
imposed mild hierarchy between mjrj and M2D- To this end, m\Y> and ni2D turn out to be 
heavier than the ones of the corresponding quarks and lie in the ranges (0.1 — 1) GeV and 
(2 - 20) GeV while the obtained M^c, M 2 $c and M 3 p are restricted to the values 10 10 GeV, 
(10 11 - 10 12 ) GeV and (10 13 - 10 15 ) GeV respectively. 
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